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Abstract 

We study the problem of finding the instability index of certain non-selfadjoint fourth order 
differential operators that appear as linearizations of coating and rimming flows, where a thin 
layer of fluid coats a horizontal rotating cylinder. The main result reduces the computation 
of the instability index to a finite-dimensional space of trigonometric polynomials. The proof 
uses Lyapunov's method to associate the differential operator with a quadratic form, whose 
maximal positive subspace has dimension equal to the instability index. The quadratic form 
is given by a solution of Lyapunov's equation, which here takes the form of a fourth order 
linear PDE in two variables. Elliptic estimates for the solution of this PDE play a key role. We 
include some numerical examples. 

1 Introduction 

The stability of steady states is a basic question about the dynamics of any partial differential 
equation that models the evolution of a physical system. Frequently, the first step is to linearize 
the system about a given equilibrium. Linearized stability is determined by the spectrum of the 
resulting differential operator A. If A has discrete spectrum, an important quantity is the instability 
index, k(A), which counts the number of eigenvalues in the right half plane (with multiplicity). 

In order to numerically evaluate the instability index of a given differential operator, its com- 
putation should be reduced to a problem of linear algebra. Particularly for problems with periodic 
boundary conditions, it seems natural to restrict A to a finite-dimensional space of trigonometric 
polynomials. Under what conditions can k(A) be computed from the resulting finite matrix? One 
difficulty is that the entries of the infinite matrix corresponding to the differential operator A grow 
with the row and column index, so that any truncation is not a small perturbation. 

If A is a self-adjoint semi-bounded differential operator of even order, then the computation of 
its instability index is well-understood through the classical work of Morse lfl8l who solved this 
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problem completely in the space of vector functions in one independent variable. The instability 
index of A agrees with the dimension of the positive cone of the corresponding quadratic form. It 
is invariant under congruence transformations that replace A with T*AT. The instability index can 
be estimated by variational methods, or computed directly from the zeroes of the corresponding 
Evans function. 

Understanding the spectrum of a non- self adjoint operator is a much harder problem. It is not at 
all obvious how to restrict the computation of its instability index to a finite-dimensional subspace, 
or how to even estimate its dimension. Furthermore, the numerical calculation of eigenvalues can 
be extremely ill-conditioned even in finite dimensions. One impressive example is the matrix 

10 6 

A= I 10" 2 





-(io 6 ) 

The Mat lab function eig(A) gives for the eigenvalues the numerical results Ai = —0.8, A 2 /3 = 
2.4 ± 1.7 i, which suggests an instability index of k(A) = 2. However, the accuracy of the compu- 
tation is poor. Denoting by V the matrix that contains the (numerically computed) eigenvectors in 
its columns, and by E the diagonal matrix that contains the (numerically computed) eigenvalues, 
then 

norm {A - VEV^) = 7.6 . 
On the other hand, A is similar to an upper triangular matrix 



A = T 2 10 6 T , where T 





and we see that actually Ai = A2 = 1, A3 = 2 and k(A) = 3. In contrast, the eigenvalues of the 
symmetric matrix 

/ 10 4 + 1 10 6 10 4 
B = 10 6 2 -(10 6 ) 

\ 10 4 -(10 6 ) -(io 4 -i) 

can be determined with the much better computational accuracy 



norm (B - V EV~ l ) = 3.6 * 10" 
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Note that B differs from A only in the signs of two off-diagonal entries. The chance of encounter- 
ing a matrix with moderately- sized entries and a badly conditioned eigenvalue problem increases 
rapidly with the dimension of the matrix (see lfi~4l 1241 ). Such examples demonstrate that the sta- 
bility problem for a non-selfadjoint operator cannot be easily solved by direct computations of the 
spectrum. 

In this paper, we examine the computation of the instability index for differential operators of 
the form 

A[h] = -h"" - (a(x)h)"+ (b(x)h)' - c(x)h , (1.1) 
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acting on 27r-periodic functions. Such operators appear as linearizations of models for thin liquid 
films moving on the surface of a horizontal rotating cylinder. The resulting flows are called coat- 
ing, if the fluid is on the outside of the cylinder, and rimming, if the fluid is on the inside of a 
hollow cylinder. They appear in many applications, including coating of fluorescent light bulbs 
when a coating solvent is placed inside a spinning glass tube, different type of moulding processes 
and paper productions. 

One would expect the flow to become unstable, if the fluid film is thick enough so that drops of 
fluid can form on the bottom of the cylinder (in case of a coating flow) or on its ceiling (in case of 
a rimming flow). In both cases, surface tension and higher rotation speeds should help to stabilize 
the fluid, but may also allow for more complicated steady states. 

The operators in Eq. dl . lb appear as linearizations of the flows about steady states, when the 
dependence on the longitudinal variable in the cylinder is neglected. Benilov, O'Brien and Sazonov 
||5l studied the convection-diffusion equation 

d /, . dh\ 



A[h] = —fh + e sin x—) 

dx\ dxJ 



with periodic boundary conditions on [0, 2%], which corresponds to a singular limit of a rimming 
flow where surface tension is neglected. This operator has remarkable properties: For \e\ < 2, all 
its eigenmodes are neutrally stable, but the Cauchy problem 

jh = A[h], h(0) = h 

is ill-posed in any Sobolev and Holder space of 2ir -periodic functions. The underlying cause is the 
sign change of the diffusion coefficient as a; — > x + ir. This phenomenon of explosive instability of 
a system with purely imaginary spectrum was studied analytically by Chugunova, Karabash and 
Pyatkov [|9]|, who explained it in terms of the absence of the Riesz basis property of the set of 
eigenfunctions. The spectral and asymptotic properties of A are of interest in operator theory and 
were analyzed in Ifl2ll26ll8l[l0l 

One should expect the explosive instability to disappear in complete models that includes the 
smoothing effect of surface tension. Such models have been proposed, for example, by rt20ll2TTl . 
In [Vj, the authors linearized this model about some approximation of a positive steady state 
solution to obtain 

. rjl/ , d f ., dh ( dh d 3 h\} 

A[h](x) = -- |(1 - ttl oosx)h + a 2 smx- + a 3 + — j j (1.2) 

with periodic boundary conditions. Here, the parameter ct\ is related to the gravitational drainage, 
a 2 is related to the hydrostatic pressure (in lubrication approximation model this coefficient is very 
small), and the parameter a 3 describes surface tension effect. They showed numerically that a 
sufficiently strong surface tension can stabilize the film provided that the other coefficients are not 
too small. For smaller values of a\ and a 2 , capillary effects destabilize the film. The number of 
unstable eigenvalues of A grows if a 3 is decreased. 
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We will consider operators given by Eq. (11.11) acting on L 2 [0, 2%] with periodic boundary condi- 
tions. We assume that the coefficients a(x), b(x) and c(x) are bounded smooth periodic functions. 
We will show that the instability index of A is determined by its projection to a sufficiently large 
finite-dimensional subspace of L 2 . The dimension of the space depends on a suitable norm of the 
distributional solution U of the partial differential equation 

A*U(x,y)=5 y . x (1.3) 

with periodic boundary conditions on [0, 2%] x [0, 2ir). Here, the differential operator A is defined 
by applying the single-variable differential operator A to the functions F(-,y) and F(x, ■) and 
adding the results; symbolically 

We note that Eq. (|1.3I) has a unique solution if the spectra of A and —A* are disjoint [2]. Let 
Uo(x, y) be the solution of Eq. (11.31) with a(x) = b(x) = and c(x) = 1. We will see below 
that U is piecewise smooth, with a jump in the third derivative across the line x = y, and that 
U(x,y)-U (x,y)eH\ 

To describe our results, denote by P/v the standard projection onto the space of trigonometric 
polynomials of order N, 

p N [<i>){x)= k*v px - (i.4) 

\p\<N 

In Proposition |7J] we show that 

k(A) = k(P n U- 1 P n ), (1.5) 

provided that 

N 2 >2M{1 + \\U(x,y) -U (x,y)\\ Hi ). 

The constant is given by 

M = 0.52 (\\a\\ H i + \\b\\ H i + \\c - l\\ H i) . (1.6) 

The significance of Eq. (11.51) is that it allows to compute the instability index of A from the finite 
matrix that describes the restriction of U to the finite-dimensional subspace 

(Range (I - P)) ±u = Nullspace (P N U) . (1 .7) 

The weakness of this result is that both the condition on iV and the computation of the subspace 
involve the unknown function U, which is defined as the solution of a partial differential equation. 
The existence of such a solution, and its norm, depend sensitively on the spectrum of A, which is 
exactly the unknown quantity we are concerned with. 

It is tempting to consider instead the matrix obtained by truncating the Fourier representation 
of A at a suitable high order N. Our main result, Proposition 17 A[ guarantees that 

k(A) = k(P n AP n ) . (1.8) 



4 



provided that 

N 2 > M(l + V / M+I) (1 + | \U(x, y) - U (x, y) | \ n .) . 

Note that only the norm of the unknown function enters into the condition on N, and that the 
identity in Eq. (11.81) does not involve U at all. 

The selection of N and the problem of estimating this norm will be discussed at the end. 

Let us add a few words about the proofs. Our analysis relies on the indefinite quadratic form 
defined by the self-adjoint operator U. Classical results, which will be discussed in the next section, 
state that 

k{A) = k{U) , 

and that the positive and negative cones of U contain the invariant subspaces associated with the 
spectrum of A in the right and left half planes, respectively. The key to Eq. (11.5b is that the quadratic 
form is negative on high Fourier modes, because the fourth order term in A dominates the lower 
order derivatives. As part of the argument, we derive an addition formula for the instability index 
of a self-adjoint operator in terms of its restriction to suitable subspace. The proof of Eq. (|1.8b 
combines Eq. (11.5b with estimates for the off-diagonal terms in the Fourier representation for U . 

One of the possible extensions of our results could be an application of a similar method to 
obtain the estimations on the size of the finite dimensional truncation in the case of a more general 
forth order differential operator with the third order derivative term which is absent in 1 1.21 

2 Lyapunov's equation 

The partial differential equation (11.3b is an instance of Lyapunov's equation 

A*U + UA = V, (2.1) 

which was first considered by Lyapunov in the case where A and U are n x n matrices, and V 
is symmetric and positive definite. (In Eq. (11.3b , V = I.) Assuming that a symmetric matrix U 
solves Eq. (12.1b , Lyapunov proved that all eigenvalues of A have negative real part, if and only if 
U is negative definite. The follwoing generalization is due to Taussky [|22l . 

Theorem 2.1 (Taussky). Let A be an n x n complex matric with characteristic roots ctj, with 
+ o>k 7^ Ofar all i, k = 1, . . . n. Then the unique solution U of Lyapunov's equation with V = I 
is nonsingular and satisfies k(U) = k(A). 

The problem of obtaining information about the sign of eigenvalues of A in situations where 
both V and U may be indefinite and have non-trivial kernels remains an area of active research. 

Lyapunov's equation has many applications in stability theory and optimal control. In typical 
applications, k(A) = 0, so that the system is asymptotically stable, and U is used to estimate the 
rate of convergence. Eq. (12.1b is a special case of Sylvester's equation 

AX — XB = C , 
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which has been studied extensively in Linear Algebra, Operator Theory, and Numerical Analysis. 
It is known to be uniquely solvable, if and only if the matrices A and B have no eigenvalues in 
common. In particular, Eq. (12.11) has a unique solution if the spectra of A and —A* are disjoint. 
Since V is self-adjoint, a unique solution U is automatically self-adjoint as well. These results were 
extended to bounded operators on infinite-dimensional Hilbert spaces by Daleckii and Krein IfTTTl 
and to unbounded operators by Belonosov (21 El. 

Before stating Belonosov's result, we recall that a closed densely defined operator A on a 
Banach space is sectorial, if the spectrum of A is contained in an open sector 

S={zEC ||arg(A -*)| < 

with vertex at A G M and opening angle 9 < tt/2, and the resolvent R\(A) — (A — A/) -1 is 
uniformly bounded for A outside S. Sectorial operators are precisely the generators of analytic 
semigroups. The sector S is invariant under similarity transformations, and does not change if the 
norm on the space is replaced by a equivalent norm. 

Theorem 2.2 (Belonosov). Let Abe a sectorial operator on a separable Hilbert space H. Assume 
that 

a{A) n <t(-A*) =0. 

Then for any bounded operator V on H, the Lyapunov equation ({2.1\) has a unique solution U in 
the class of bounded operators on H. Then U is invertible in the general sense, i.e. its inverse is 
densely defined but can be unbounded operator 

k(A) = k(U) = nil/- 1 ). 

Belonosov actually proved more general existence and uniqueness results for the Sylvester's 
equation in Banach spaces. 

To explain the geometric meaning of Lyapunov 's equation, we introduce on H the indefinite 
inner product 

[0,V] = W,V>>- (2-2) 
If U has trivial nullspace and k(U) < oo, then H equipped with [-, •] is called a Pontryagin space, 
and will be denoted by II. The concepts of orthogonality and adjointness are defined on II in the 
natural way with respect to the indefinite inner product [x, y] . A subspace X C II is called positive 
if [/, /] > for for every non-zero vector f E X, and negative if [/, /] < for every non-zero 
f E X. Maximal positive subspaces have dimension k(U), while maximal negative subspaces 
have codimension k(U). 

Let 4>(t) = e tA (f)o be a solution of the evolution equation 

j t <i>(t) = A<f>(t) , 0(0) = 0o- 

Lyapunov's equation guarantees that the value of the quadratic form Q{4>) = [0, 0] strictly in- 
creases with t, 

iQW)) = {{A*U + UA*)<P(t), 0(t)> = <V0(t), 0(t)> > . 
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Denote by M + (A) the invariant subspace associated with the part of the spectrum of A located in 
the right half plane. If G M+(A), then 

Q(<t>) > lim Q(e tA cf) =0, 

which shows that M + (A) is a positive subspace of IT. A similar argument with t — ► oo shows 
that the complementary subspace M_(A), which corresponds to the spectrum of A in the left half 
plane, is a negative subspace of IT. Since Eq. (12.11 ) excludes purely imaginary eigenvalues, these 
subspaces are maximal, and consequently k(A) = k{U). 

One can also interpret Lyapunov's equation as a dissipativity condition on A with respect to 
the Pontryagin space II. In general, a densely defined linear operator A on II called dissipative if 
Re [Af, f] < for all / G Dom(^4). It is maximally dissipative if it has no proper dissipative 
extension in II. Assuming Lyapunov's equation, we compute for <p ^ 

Re [Acf>, 4>] = Re (UA<f>, 4>) = 1 -{(A*U + UA)<f>, 4>) = ^(Vc/>, 0) > , 

i.e., —A is dissipative. In this framework, the analogue of Belonosov's theorem was proven by 
Azizov [QQl (but note that Azizov formulates the result in terms of Im rather than Re ): 

Theorem 2.3 (Azizov). Let A be an operator on IT such that —A is maximally dissipative. Then 
there exist a maximal nonnegative subspace U + and a maximal nonpositive subspace Il_ of TI such 
that 

Rea(A| n+ ) > 0, Re a(A\ n _) < . 
Moreover, we can choose II + and Il_ to be invariant subspaces for A, and 

n + d m + {A) , n_ d m_(A) . 

If additionally, Re [Af, f] > for all nonzero f G Dom(A), then M + (A) and M_(A) are 
themselves maximal positive and negative subspaces for IL respectively, and 

M + (A)+M_(A) = II. 

The second part of Azizov's theorem implies that k(A) = k(U) provided that V in Eq. (|2.1I) 
is positive definite. This agrees with the conclusion of Theorem 12.21 but note the difference in 
the hypotheses: Belonosov's assumption that A is sectorial provides resolvent estimates that allow 
to represent U as a contour integral (thereby proving existence), and the analytic semigroup e tA 
appears in the proof that k(A) = k(U), as sketched above. In contrast, Azizov's theorem does 
not require A to be sectorial, but starts instead from a given solution to Eq. (12.11) . In the special 
case where k(U) = 0, Theorem 12.31 reduces to a theorem of Phillips that characterizes maximal 
dissipative operators as generators of strongly continuous contraction semigroups. In particular, 
the spectrum of A lies in the closed left half plane (see 11251 , Corollary 1 in Section IX.4). 

In the case where A is a sectorial differential operator of even order on an interval [a, b] Be- 
lonosov proved that the solution of Lyapunov's equation with V = I is given by a self-adjoint 



7 



bounded operator [4J. His results are formulated for "split" boundary conditions that do not cou- 
ple the values at the two ends of the interval. Belonosov's results were extended to second-order 
sectorial differential operators with non- split boundary conditions by Tersenov [|23l . The operators 
we consider here are of fourth order with periodic boundary conditions. 

It is an interesting open question how to take advantage of the freedom to choose an arbitrary 
positive definite self-adjoint bounded operator V for the right hand side of Eq. (I2.ll ). For instance, 
if A is a sectorial non- self adjoint differential operator, can V be chosen in such a way that the 
solution U is the inverse of a differential operator? 



3 Spaces and norms 

We start with some estimates for the differential operator in Eq. (Il.lt . We will work in L 2 = 
L 2 ([0, 27r]), and will use periodic boundary conditions throughout. The inner product and norm 
are denoted by 

(f,g) = j** f{x)g{x)dx, \\fU=^ \f{x)\ 2 dx ' 

For the Fourier coefficients we use the conventions 

In the hope of minimizing confusion, we will denote functions on [0, 27r] by lowercase letters such 
as (/, <fi, . . . ), and functions on the square [0, 2tc] x [0, 27r] by uppercase letters (F, $, ....). Abusing 
notation, we will identify a function F(x, y) G C 2 with the corresponding integral operator F on 
L 2 . By Schwarz' inequality, 

\{F<P^)\= 1^ r F(x,y)<P(y)i>{x)dydx <\\F(x,y)\\ c *\\<l>\\ti>\\il>\\ti>, 
Jo Jo 

and consequently 

\\F\\ L ^ L 2<\\F(x,y)\\ C 2. (3.1) 

Operators on functions of two variables will be denoted by calligraphic letters (J 7 , Q, . . . ). Given 
a single-variable operator F, we denote by T x or T y the operators that acts on the x- or ^/-variable 
of a function y) while keeping the other one fixed. 

On the Sobolev spaces H s = H s ([0, 27r]) with periodic boundary conditions, we use the norms 



2 HS = 2nj2^+p 4 r /2 \hp)\ 2 - 



p=—oo 
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The corresponding Sobolev spaces of doubly periodic functions on [0, 2tt] x [0, 27r] will be denoted 
by Ti s , and their norms are defined by 
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4tt 2 J2 (2+p 4 + g 4 r /2 |F(p,g)| 2 . 

p,q= — oo 

Note that for s = 0, this agrees with the definition of the £ 2 -norm as the square integral. The 
choice of the Fourier multipliers (1 + p 4 ) s / 2 and (2 + p 4 + q 4 ) s l 2 in place of the standard (1 + p 2 ) s 
and (1 + p 2 + q 2 ) s allows for an easier comparison between functions of one and two variables. 
Finally, we denote by D the unique positive definite self-adjoint operator on L 2 such that 

£>V = <P"" + <f> • (3-2) 

This is a first-order pseudodifferential operator that provides an isometry from H s+1 onto H s for 
every value of s. 

The domain of the operator A in Eq. (11.21) consists of periodic functions in if 4 [0, 2tt], and its 
adjoint is given by 

A*[f} = -f""-a(x)f"-b(x)f'-c(x)f. 
In particular, A is self-adjoint, if b(x) = a'(x). 

Lemma 3.1. For any a e H 1 and every (f> 6 L 2 , we have 

\\a<j)\\ L 2 <O.52||a||^i||0|| L 2. 

In particular, 

\\A* + D*\\h2^v <M, 
where M is the constant from Eq. &1.6\) . 

Proof. Since 

oo 

£(l+p*)- 1/2 «1.68, 

p=— oo 

we have, for a £ H 1 , 

I 1 \ x l 2 
sup ||a0||x2 < ||a||ioo < ( — y (1+p 4 )^ 1 ' 2 ) | \a\ \ H i < 0.52 | \a\ \ H i . 

For the second claim, we use that for <\> E H 2 

\\(A* + D 4 )^ < ||a(x)0"|U 2 + ||^)0'llL 2 + ||(c(x)-l)0|| i2 

< ||o|| £ «,||0|| fla + ||&|U=o|M| H i + \\c- l\\ L oo\\<f>\\v 

< M\\<j>\\H*. 

□ 
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Lemma 3.2. A is sectorial. 

Proof. It suffices to show that the Hausdorff set { (/, Af) \ f e Dom(A), 1 1/| | = 1} is contained 
in a closed sector 

S={\ }U{zeC : |axg(A -z)| < 9} 

with some vertex A and opening angle 9 < n, and that A — X I is invertible (see p. 280 of IfTTl ). 
Choose 

Ao = - (l + max{— a"(x) + b'(x) — c(x)} + (max[a(x)] + ) 2 ) 

2 X X 

9 = tan (max \a'(x) — b(x)\) . 

X 

We estimate, for f E H 2 

Re (/, (Ao - A)f) = ^ \ff - a(x)\ff + (a + \a\x) - \b'{x) + c(x)) \f\ 2 dx 

^ 2 - E -^f\m\ 2 (3-3) 

p=— oo 
1.....2 



2 i m i l-H" 2 



This shows that the spectrum of A lies in the half plane Re z < \q — \. Similarly, 

|Im(/,(A -A)/)| < r \a\x)-b(x)\\lmf'f\dx 

Jo 



< 27rmax \a'(x) — b(x)\ > |p||/G°)| 2 - 

X ' J 

p=— oo 

For 1 1/ 1 1 = 1 it follows that 

|Im (/,(A -A)/)| , v / 2\p\\ 

, ? „ 77 / \ tttt - < (maxo i -6x J sup— -) = max o (x) - b(x) \ , 

Re (/, (A - A) J) x v peZ l+p 4 / x 

which yields the claim. □ 



The lemma implies that the Cauchy problem for A has a unique solution for every initial value 
ho G L 2 . This solution is analytic in t for £ > 0, and for any fixed t > 0, the function ft,(t, ■) e 
Dom(y4). If the coefficients of A are analytic, then h is analytic in both variables for t > 0. An 
application of the Lax-Milgram theorem similar to Lemma |4~T| below shows that (A — A)^ 1 maps 
L 2 into H 2 . It follows that the resolvent is a compact operator of the Hilbert-Schmidt type, and 
that the spectrum of A is discrete. 
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4 The integral kernel U (x, y) 

Let A be the differential operator from Eq. (11.11) . Theorem 12.21 implies that Lyapunov's equation 
has a unique solution U, provided that the spectra of A and A* are disjoint. Our goal is to show 
that U admits an integral representation 

U{f)= rU(x,y)f(y)dy, 
Jo 

and to derive bounds on U (x, y). Equation (12.11) requires that ((A*U + UA*)4>, ip) = (4>, ip) for all 
smooth periodic test functions functions 0, tp. This means that U{x, y) is a distributional solution 
of the partial differential equation (11.31) . 

Let us solve Eq. (11.3b in the special case A = —D 4 , given by 

A [f] = -f""-f. 

By our choice of norms, — A defines an isometry from H 4 onto L 2 . Since A has constant coeffi- 
cients, the unique solution can be written as U (x, y) = u (x — y), where 

2A Q u = 6 , 

in other words, 2U (x, y) is the Green's function of A on [0, 2ix] with periodic boundary condi- 
tions. One can compute u (x) explicitly as a linear combination 

X — IX , X — IX „ . X — 7X . X — 7X 

u Q (x) = Ci cos — =- cosh — = — h Co sin — — smh — =- , 

where the coefficients are adjusted so that u is periodic and twice differentiable, and its third 
derivative jumps by —1/2 at x = 0. From this representation, it is clear that U is smooth away 
from the line x = y, and that U G C 2 ' 1 C H 3 . Alternately, we easily obtain from the Fourier 
representation of A that Uo(p) = —^j^prj, an d 

1 °° 1 

U (x,y) = V e ip{x ~ y) . (4.1) 

ov ,yj Atx ^ 1+p 4 

p=—oo 

In particular, U (x, y) G H s for all s < ~, and | \U \ \h 3 < 1- 
It remains to analyze the difference 

K(x, y) := U (x, y) - U Q (x, y) . 

By definition, K solves the partial differential equation 

A*K(x, y) = - {A* - A ) U (x, y) . (4.2) 

The second order differential operator A* — A® maps maps H 2 into L 2 , see Lemma [3TT1 A weak 
solution of this equation is provided by the next lemma. 
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Lemma 4.1 (Construction of K). The resolvent of A* is compact and maps C 2 into H 2 . 

Proof. Let Ao be the vertex of the sector computed in Lemma [3^21 and assume that F(x, y) E C 2 . 
We verify that the equation 

(2\ -A*)K(x,y) = F(x,y) 

satisfies the assumptions of the Lax-Milgram theorem, as stated in [Evans, PDE, p. 297] lfT3ll . 
Define a bilinear form on on smooth doubly periodic functions \l/ 

J B($,^) = (($,(2A -X)^) £2 . 

Then B is extended continuously to H 2 by 

£?($, *) = ($, V) H 2 + 2A ($, - ($, (A* - Ao)*) c2 . 

On the other hand, it follows from Eq. (13.31) that 

1 



-||$ 



|2 
\H 2 ■ 



Finally, the map 

defines a continuous linear form on H 2 . The Lax-Milgram theorem asserts that there exists a 
unique function K(x, y) e H 2 such that 

B(if l «r) = (F,*) jt? 

for all ^ G 7i 2 . By the resolvent identity, the equation 

(A*-\)K(x,y) = F(x,y) 

has a unique weak solution in Ti, 2 for every value of A that is not an eigenvalue of A* and every 
F(x,y)eC 2 . □ 

Lemma 4.2. IfK(x, y) G 7i 2 so/ves £g. f[4.2|) . y) G 7i 4 , and 

||^(z,y)|| W 4 < 2M\\U (x,y) + X(x )2 /)|| w2 , 
where the constant is given by Eq. rti.(5D . 
Proof. If K(x, y) solves Eq. (14.21 ), then 

AoK = -(A*-Ao)(U + K), 

and we conclude that 

\\K(x,y)\\ H 4 < \\A*-A \\n^\\U + K\\ H 2<2\\A*-A \\ H 2^ L2 \\U + K\\ n ,. 
The proof is completed with Lemma I3TT1 □ 
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5 Estimates for the operator U 



In this section, we derive bounds for U = Uq + K as an operator on L 2 . Since K(x, y) E 7f 4 , while 
Uq(x, y) £ H s only for s < 7/2, the Fourier coefficients of K(x, y) decay more quickly than the 
Fourier coefficients of U (x, y). This in turn implies that the restriction of U to high Fourier modes 
is dominated by Uq. In this section, we provide the relevant estimates. 

As a consequence of the regularity result in Lemma H~2l we see that U defines a bounded linear 
operator from L 2 to ff 4 , with 

\\U\\ L 2^ H * < \\U \\ L 2^ H 4 + \\K\\ L 2^ H 4 <- + \\K{x,y)\\ H *. 

We have used that D A X U = \8 and applied Eq. (f3j} to D 4 K(x, y). 

One attractive property of the 7f 4 -norm is that it depends only on the magnitude of the Fourier 
coefficients, not on the phases. In contrast, the operator norm 

oo 

"-'| M 4= sup (£> 4 F0,V)=47r 2 (l+p 4 )Hp,q)i(q)j>(p) 



\\<P\\ L 2=\W\ L 2-- 



p,q=—oo 



can change drastically if we replace F(p, q) by \F(p, q)\. This dependence on cancelations can 
cause difficulties in estimates: Multiplying the Fourier coefficients of F with factors a(p,q) E 
[0, 1] will not necessarily decrease the operator norm. On the other hand, the if 4 -norm provides 
only a rather loose bound on the norm of the corresponding integral operator. For instance, the 
kernles U (x, y) (and consequently U (x, y) does not lie in H 4 , even though | |f7 | \tf-+H* — §• 

We find it useful to introduce another norm on integral kernels that lies between the H 4 -norm 
(as a function of two variables), and the operator norm (as a linear transformation from L 2 to if 4 ). 
By construction, this norm depends only on the modulus of the Fourier coefficients. 

Lemma 5.1 (Auxiliary norm). Define, for smooth doubly periodic functions F 

oo 

|||F|||:=47r 2 sup V (2 + p 4 + q A )\F(p, q )\\4>(p)\\^( q )\. 
\\<t>\\=U\\=i, ^ 



p,q= — oo 



Then 
and 



\F\\\<\\F(x,y)\\ n ,, 



\\\F\\ \ > max{||F|| L 2_ H 4, \\F\\ H -4_, L 2, 2||F||. ff -2_ >iJ -a} . 
Proof. From the Fourier representation, we see that 



oo oo 

< sup _ 4vr 2 E ^ + v A + q A )\Hv,q)\\Hv,q)\ 

\\^{ 3S >V)\\c i ~ 1 p,q=ioo p,g=— oo 



< sup (AqF, $) £ 2 

ll < f(^,S/)ll £ 2=l 

= \\F(x,y)\\ H *. 



13 



On the other hand, 



oo 



\\F\\v>->h* = \ \D 4 F\\ L 2^ L 2 = sup V (l+p 4 )0W(g)F(p,g) < 
and similarly 

\\F \\h~ 2 — >H 2 — ~\\\F\\\ ; I \F\ \h~ 4 — *l 2 — 1 11-^1 1 1 • 



□ 



We note that if F has positive Fourier coefficients, then \ \\F\\ \ agrees with the operator norm 
of AqF as a linear transformation from L 2 into itself. In particular, 1 1 \Uo\ \ \ = 1. 

Lemma 5.2 (Tail estimate). Assume that K(x, y) solves Eq. {\4.2\) , and let M be given by Eq. {\1.6\) . 
Then 

\\\K - P N KP N \\ \ < MN- 2 (\\\U + K\ 
Proof. Using Eq. (14.21) together with the definition of the norm, we obtain 

\\\(K - P N KP N \\\ = Arr 2 sup 1^)11^)11(^^)^9)1 

11011 = 11^11=1 | p |>7V or |g|>7V 



;i +p 4)l/2 + ( 1+9 4 V l 2 

\p\>N or |g|>iV 

< mn-' 2 (\\\u + k 



□ 



Lemma 5.3. If U solves Eq. ([O]), then \ \\U\\ \ > I. 
Proof. Write U — U + K, and estimate 

\\\U\\\>\\\(I-P N )U (I-P N )\\\-\\\(I-P N )K(I-P N )\\\. 

The first summand is bounded below by 1 because AoUo = I, and the second summand is bounded 
by MiV~ 2 |||[/||| according to Lemma E21 We conclude that (1 - MiV- 2 ) |||[/||| > 1 for each N, 
and the claim follows upon taking N — > oo. □ 



6 Addition rule for the instability index 

We return to the Pontryagin space II introduced in Section [2J with the indefinite inner product 
given by Eq. (12.21) . Let LIi be a finite-dimensional subspace of II, and let 



IL 



11^ = {/ell [f,g] =0 for all g eU,} 
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be its [/-orthogonal complement. By construction, dimlli = codimn 2 . The natural question is 
can we compute k(U) from the restrictions fv(£7|rii) and k(U\u 2 )^ The difficulty is that n need not 
be a direct sum of 111 and n 2 , because the two subspaces may intersect non-trivially in a subspace 
where the quadratic form vanishes. 

A subspace X c II is called neutral, if [0, 0] = for all (ft £ X. Two finite-dimensional 
neutral subspaces X and Y of H are ITskewly linked, if 

dimX = dimF 

and the inner product [., .] does not degenerate on the direct sum X+Y. In particular, no vector of 
X different from is orthogonal to the skewly linked subspace Y, and vice versa. 

Theorem 6.1 (Theorem 3.4 L16J). Let 111 be an arbitrary subspace ofH, let II2 be its U -orthogonal 
complement, and let X = Hi D il 2 be their intersection. There exists a neutral subspace Y C II 
that is skewly linked to X and provides a U -orthogonal decomposition 

n = n; © (x+y) © n' 2 (6.i) 

where 

n x = n; © x, n 2 = n 2 © x. 

The theorem was originally formulated for the case of regular Pontryagin spaces, where the 
quadratic form U is a bounded operator with bounded inverse. Under the assumption that 111 is 
finite-dimensional, the result easily extends to the situation where the inverse of U is unbounded 
but densely defined. Although the above decomposition is not unique in general, it yields the 
following addition formula for instability indices: 

Proposition 6.2. Let Hi be a finite-dimensional subspace of II, and let H 2 be its U -orthogonal 
complement, Then its instability index is given by 

K (U) = k{U\ Ui ) + K(U\n 2 ) + dim(IIi n II 2 ) . 
Proof. Theorem 16. ll provides subspaces H[ and W 2 such that 

k{U) = k{U\ UIi ) + k(U\w 2 ) + k(U\ x+y) ) . 

By construction, we have k(U\ji[) = «(Z7|ni) and k(U\u' 2 ) = k(U\u 2 )- It remains to compute 

k(U\ x +y))- 

Since X and Y are skewly linked and finite-dimensional, there exists for each basis (pi, 2 , ...(j) m 
of X a basis ipx, ip 2 , •••'0m of Y such that [0j, ipj] = 5ij = 1, m). By expanding an arbitrary 
element h 6 X+Y as 

m m 

i=i 3=1 
the indefinite inner product can be expressed as 

m 1 / m m 

[h, h\ = 2 «iA = 2 Y,^ + - - Prf 

i=l \i=l i=l 

This is an explicit representation of the indefinite inner product in terms of positive and negative 
squares, which shows that k(U\ x +y) = dim(X). □ 
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7 Restriction to finite dimensions 

We first prove the claim in Eq. (11.51) . 

Proposition 7.1 (Projecting out high Fourier modes). Let A be given by Eq. ([777]). Assume that 
the spectra of A and —A* are disjoint, and let U (x, y) be the kernel of the unique solution of 
Lyapunov's equation was constructed in Section^ If 

iV 2 >M|||f/|||, (7.1) 

where M is the constant from Eq. ( |i.<5l) . then 



k(A) = k(P n U- l P 



N) 



Proof. By Theorem 12.21 we have k(A) = k{U). Let = (U(j),ip) be the indefinite inner 

product associated with U. Choose n 2 to be the range of I — P N , and let Hi = n 2 Lc/ be its 
U -orthogonal complement. We will show that 

k(U) = k(U\ Ui ) ■ (7.2) 
This will establish the conclusion, because 

and U maps tlx isomorphically onto the range of P N . 

Let us apply [-, ■] to D 2 (f), where <fi e H 2 and D is given by Eq. (13.21) . Writing U = U + K, 
and using that D 2 U D 2 = — |j, we see that 

[D 2 ^D 2 <P} = ((U + K)D 2 <P,D 2 <P) = — ^||0|| 2 + (D 2 KD 2 (f), <p) . 
We replace 4> with (I — Pn)4>, and use Lemma 15721 to obtain 

[D\l-P N )4>,D\l-P N )<f>] < ~(1-\\\(I-P N )K(I-P N ) 

< ~(l-e N )\\(I-P N )^ 3 
where e N = MN~ 2 1 1 \U\ 1 1 < 1. It follows that 



U\ n2 <(l-s N )U o \ n2 <0 (7.3) 

as quadratic forms on H 2 . In particular, k(U\u 2 ) = 0, IIi fl n 2 = 0, and Eq. (|7.2I) follows with 
Theorem 16. II □ 

For our final result, we want to replace TIi by the range of the projection P^ from Eq. (11.41) . 
The next two lemmas concern the restriction of U to the range of P^. 
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( A n 


A 12 \ 


■ ^=1 


f Un 


C/12 \ 


V A 21 


A 22 J 


V ^21 


U22 ) 



Lemma 7.2 (Lyapunov equation for PnAPn). Let A be given by Eq. f li.il) . Assume that the spec- 
trum of A and —A* are disjoint, and let U (x, y) be the kernel of the unique solution of Lyapunov' s 
equation that we constructed in Section\4\ If 

iV 4 > M 2 \\\U\\\, 

then 

{P N AP N )*{P N UP N ) + (P N UP N )(P N AP N ) > c N P N , (7.4) 
where cn = 1 — M 2 N~ A \ \ \U\ \ \. In particular, 

k(P n AP n ) = k(P n UP n ) . 

Proof. For G L 2 , we write <pi = Pn4>, 4>2 = (I — Pn)4> an d decompose 

From Eq. (12.11) . we see that Un solves Lyapunov's equation 

A* u U n + U n A n = V 

with V = In — A\ 2 U 2 i — U12A21. We claim that the right hand side is positive definite on the 
range of P N . 

To prove this claim, first observe that we can replace A by A — A and U by K = U — U in 
the definition of V, because A and Uo are diagonal in the Fourier representation. We estimate 

\\A* 12 U 21 + Ui 2 A2i\\ L ^ L 2 = \\P N (A* - A )(I - P N )KP N + P N K(I - P N )(A - A )P N \\ L ^ 

< sup { P > 4 \ T M\\\K-P N KP N \\\ 

\p\>N or \q\>N ±+P + q 

< M 2 N~ 4 

by Lemma 15721 It follows that V > cnPn > as quadratic forms on the range of P N . Since 
PnUPn is a finite matrix, the conclusion of the lemma follows with Taussky's theorem. □ 

Lemma 7.3. Under the assumptions of the previous lemma, PnUPn is invertible on the range of 
Pn, and 

\\(D 2 P N UP N D 2 y 1 \\ L ^ L ,<2 1 -±^. 

c N 

Proof Let us write A N = D- 2 P N AP N D 2 and U N = D 2 P N UP N D 2 . By Eq. £774]), we have 

A* N U N + U N A N > c N D A P N 

as quadratic forms on the range of P N . Here, = 1 — M 2 N~ 4 \ \ \U\ 1 1, as in Lemma [7721 We apply 
this inequality to an eigenfunction O of Un 

2a n2. 



2A Re {A* N <f) , O > > c N (D 2 </> , D 2 ^) 
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where Ao is the corresponding eigenvalue. Writing O = -D 2 ip , and using once more Lemma [5^21 
we conclude that 

1 1 c/jv 1 1 i,2->z2 < ^ sup I, )J|2 < — H^- L> ||l 2 ->l 2 < L- 



^ eL 2 Ctv| |^||22 Cat Cat 

□ 

We are finally ready for our main result. 

Proposition 7.4 (Projection onto trigonometric polynomials). Let Abe a differential operator 
given by Eq. di.il) . Assume that the spectra of A and —A* are disjoint, and let U (x, y) be the 
unique weak solution ofEq. &1.3\) in 7i 2 . If 



N 2 > M(l + VT+M)|||?7||| , (7.5) 
where M is given by Eq. rti.<5D . then 

k{A) = k(P n AP n ) . 

Proof. Since U solves Lyapunov's equation, Theorem 12.21 implies that k(A) = k(U), and we 
already know from Lemma [731 that k(P n AP n ) = k{P n UP n ). We want to apply Proposition [6]2] 
in the case where IT.! is the range of P N . 

Since \ \\U\\ \ > 1 by Lemma l53l our assumption implies that 5n = MN~ 2 < 1. On 

n 2 = Rmge(P N ) ±u = {4>eL 2 \ Unfa + U 12 <f> 2 = 0} , 

we compute for the indefinite quadratic form 

[(f), 0] = (Unfa, 0l) + (Unfa, <h) + (^2101, <h) + (^2202, 02> 
= -(^21^n 1 f^l202, 02) + (^22 02, 02) • 

By Eq. (17.31) of Proposition 17.11 the last term is negative on the nullspace of P N , and satisfies the 
bound 

g ^<- i-Mira (J -P„) 

as quadratic forms. To estimate the other summand, Lemma I5T21 yields 

||^21 J D 2 ||L^<^|||(/-P)t/P|||<y 



and analogously 

WD 2 U l2 D 2 \\ L ^ L ,<^ 



The middle factor is controlled with Lemma 1731 by 



\(D 2 P N UP N D 2 y 1 \\ L ^ L ,<2- ' U 



K\\\u\W 
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We arrive at 

D*{-U 21 U^U 12 + U 22 }D> < -1(1 - 6 N \\\U\\\ - ! ^1^, (1 + M) . 

as quadratic forms. Since <%-| | \U\ | | 2 (1 + M) < (l — 5 N \ \ \U\ | |) 2 by assumption, U is negative 
definite on U 2 . It follows from Proposition l6.2l that k(U) = k(PnUPn), completing the proof. □ 



8 Numerical examples 

Before we look at examples, a few words about how to verify the hypothesis on N in Eq. (|7.1I) 
or Eq. (17.51) . The conditions involve the solution of the partial differential equation 1 1.31 A useful 
consequence of Lemmas [5T2l and l5.3l is that for 5 m := MN 2 < 1, 

1<III^III< ,-^-(1 + 111^111) ■ (8-D 
This follows by using the triangle inequality 

\\\U + K\\\ < 1 + \\\P N KP N \\\ + \\\K-P N KP N \\\ 

and solving for P^KPn in Lemma I5T21 

We propose two ways to estimate the size of 1 1 \PnUPn\ \ \. 

• Solve the partial differential equation (11.31) by a Galerkin approximation, and use this solu- 
tion to compute, approximately, the value of 1 1 \U\ \ \. If Eq. (17.11) is satisfied for some value of 
N much below the dimension of the Galerkin approximation, we can apply Proposition 17. 1[ 
and restrict U to the subspace in Eq. (|1.7I) . A basis for this subspace can be computed by 
using the Gram-Schmidt algorithm, with the inner product replaced by the indefinite inner 
product associated with U . If even Eq. (17.51) can be satisfied, then we can just restrict U to 
the range of P N . 

• Start with a value of N such that 5n = MN~ 2 < 1. Write the matrix P^APn in the Fourier 
representation, find its eigenvalues, and bring it into triangular form. Solve Lyapunov's 
equation 

P N A*P N U N + U N P N AP N = I 

for Un. In the Fourier representation, Un is a finite matrix. Compute \ max , the largest 
eigenvalue of the matrix 



{(2 + p* + q*)\U N (p,q)-l\) 



H\g\<N ■ 



Then (1 + A maa; )/(1 — 5^) is our best estimate for |||{7|||. If the condition in Eq. 17.51 holds 
with the current value of N, we are satisfied and accept the value of n(A N ) as the instability 
index for A. Else, we increase N accordingly, and repeat the above steps. 
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K(L) 




Figure 1: On the left the parameters: a\ = 0.0, a 2 = 1 and a 3 = 0.02, resulting in k ~ 190; on 
the right: oi\ = 0.0, a 2 = 1 and « 3 = 0.002, resulting in A; ps 1875. The dashed line is showing 
the suggested cut-off for the dimension of the finite dimensional subspace which is based on the 
above estimations. 



Proposition 17.41 reduces the computation of the stability index of A to a finite-dimensional 
linear algebra problem. This is illustrated in Fig. 1 for the particular example the operator A from 
0, see Eq. (11.21) . In this example, we have 



OL 2 

a(x) = 1 H sinx, b(x) 

a 3 



1 — («i + a 2 ) cosx 

«3 



c(x) = . 



In place of the constant in Eq. (|1.6I) we use the slightly smaller value 

1 + ax + 2a 2 



M= {Hp)\ + \Kp\ + \c(p)\)=2 + 



p=—oo 



«3 



The results of our computations are shown in Figure 1 . We see that if the parameter a 3 is small, 
then the surface tension is not strong enough to overcome the gravity and the model is unstable 
with the number of the unstable eigenvalues growing as the parameter a 3 decreases. 
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